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Abstract 

We propose an accurate and efficient method to compute vibra- 
tional spectra of molecules, based on exact diagonalization of an al- 
gebraically calculated matrix based on powers of Morse coordinate. 
The present work focuses on the ID potential of diatomic molecules: 
as typical examples, we apply this method to the standard Lennard- 
Jones oscillator, and to the ab-initio potential of the H2 molecule. 
Global cm~^ accuracy is exhibited through the H2 spectrum, obtained 
through the diagonalization of a 30 x 30 matrix. This theory is at the 
root of a new method to obtain globally accurate vibrational spectral 
data in the context of the multi-dimensional potential of polyatomic 
molecules, at an affordable computational cost. 

Keywords: vibrational spectra, algebraic method, Morse oscillator, anharmonic- 
ity. 

1 Introduction 

Extremely accurate rovibrational levels are obtained routinely for polyatomic 
molecules by modern spectroscopical techniques. Despite intense effort and 
significant progress, the accuracy and reliability of theoretical calculations 
is still lagging far behind expecially in the region of high overtones. Many 
present-day quantum-chemical calculations of spectroscopical properties are 
based on some power expansion of the adiabatic potential energy surface 
(PES) at its minimum [U El [3], H]. Accurate treatment of the kinetic and 
second-order potential term provides the harmonic normal coordinates and 
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frequencies. Higher-order (cubic, quartic) power corrections usually account 
for anharmonic effects. The "force field" parameters are often derived by ab- 
initio methods, or adjusted to experimentally observed levels. Either way, 
simple power terms in the normal coordinates fit very poorly the global shape 
of the actual adiabatic potential energy surface (PES) away from its mini- 
mum [S]. As a result, methods based on such power expansions retain spec- 
troscopical accuracy only up to few thousand cm~^ [HI El El El in ordinary 
simple molecules such as H2O, or H2C2; the accuracy rapidly deteriorates in 
the spectral region where dissociation is approached. 

To improve this situation, expansions of the PES in suitable internal coor- 
dinates ^lOj, for example of the Morse type [HI 1121 ES] , have been proposed. 
In this class of schemes, difficulties arise due to the generation and the need 
to diagonalize very non-sparse matrixes, which impose uncontrolled approx- 
imations, involving the neglect of a huge number of matrix elements [11]. To 
improve this state of affairs here we propose a method, based on the Morse 
oscillator, which combines the accuracy of a theory which includes anhar- 
monicity right from the start with the simplicity of the harmonic expansion 
(basis completeness, simple algebraic rules for the matrix elements, sparse 
coupling matrixes). The present work focuses on the one-dimensional (ID) 
potential, representing diatomic molecules. In this context where trivial ex- 
act numerical solution of the vibrational problem is available, the theory is 
formulated in its full detail, and its accuracy is tested against the PES of re- 
alistic diatomic potentials. This theory is constructed with careful attention 
to make it suitable as the basic building block for the successive extension to 
the multi-dimensional potentials of polyatomic molecules, which will be the 
subject of future work. 

Section [2] introduces the formalism, whose application is then investigated 
in Sec. [SI Conclusions are then drawn in Sec. HI 



The Morse potential was introduced [T3] with the purpose of producing an 
accurate description of typical spectra of diatomic molecules through an an- 
alytically solvable Schrodinger problem. Through standard manipulations of 
the three dimensional (3D) spherical problem, we come to what is usually 
named Morse equation 
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where /i is the reduced molecular mass, and the Morse Potential function 



Vm{x) = Vo{[v{x)]^ -1} , with = e-"("-"o) - 1. (2) 

This potential function is parameterized by the three quantities Vq (depth of 
the well), a (shape parameter), and Xq (position of the minimum). Vm reflects 
the main features of a diatomic molecular potential, namely one well defined 
minimum, a finite dissociation energy for large x, a strongly repulsive region 
at a; <C Xo, a finite number L of bound states and a continuum of scattering 
states. The Morse potential represents an intrinsically anharmonic oscillator, 
and indeed it was introduced originally to deal with anharmonicity exactly: 
it has the great advantage that the corresponding Schrodinger equation ([T]) 
is solvable exactly in ID, and approximately when it represents the radial 
equation of a 2-3D problem. 

In terms of the dimensionless parameter 
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the bound-state eigenfunctions can be written as 

*„(x) = iV„e-^(^)/2y(x)^-"L2(-")(y(x)) , (4) 

where 

l/(x) = (2s + l)e-"("-"«). (5) 
^ni^) ~ ^"nf ^(-^^^"^^) is a generalized Laguerre polynomial, and 



T{2s-n + l) 
2{s -n)T{n + l) 



(6) 



is a normalization constant. For integer n = 0, 1, 2, . . . [s] {[s] indicates the 
largest integer < s), the square- integrable wavefunctions "ifnix) represent the 
[s] + 1 bound states, for the specific physical parameters of the system. The 
corresponding energy eigenvalues are 

W^{x) = -^{s-nr. (7) 

In addition to its many advantages, the Morse potential has a few draw- 
backs: mainly, it does not approach the dissociation limit (here E = 0) 
with a power law as realistically expected, and it does not quite diverge at 
x = 0. On the other hand, a more realistically accurate potential V{x) can 
be obtained by adding a suitable correction Vd{x) = V{x) — Vm{x) to the 
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starting Morse potential. This correction is often written as an expansion, 
whose detailed form is determined according to some convenient recipe. Sev- 
eral approaches in this direction are documented in the literature, from the 
perturbed Morse oscillator model (PMO) [lEl [171 [18], to models based on 
algebraic methods [IHl [20l [21] . The PMO approach is based on an expansion 
of Vd{x) in a power series of v{x) 

Vd{x)=J2(^^ H^)f^ (8) 
i=4 

then treated perturbatively. In contrast, algebraic approaches based on dy- 
namical symmetries lead to an expansion of the Casimir operator of a Lie 
group, representing the whole Morse Hamiltonian [191 EDI [211 122]. Unfortu- 
nately, as this operator is of mixed kinetic and potential character, it is hard 
to relate this second type of expansion to a well defined molecular problem. 

Here we propose a substantial improvement of the PMO. This method is 
based on an expansion of the form 

^ max 

Va{x) =J2a, [v{x)]\ (9) 

i=3 

including powers of v{x) up to order Nmax, usually a small integer, e.g. 4 or 10. 
The differences with the PMO are two: (i) the inclusion of the [f (x)]^ term, 
which the PMO approach must omit due to the need to limit the number of 
low-order terms contributing to that perturbative scheme, and which allows 
much better flexibility and accuracy, as shown is Sec. [31 (ii) the exact, rather 
than perturbative solution of the resulting quantum mechanical problem. To 
this respect, the implementation of algebraic rules for the calculation of the 
matrix elements allows us to treat the third-order term on the same footing 
as all other terms in V^: we compute their matrix elements, and subsequently 
diagonalize a sparse matrix. 

The choice ([9]) of the form of the potential correction realizes a satisfac- 
tory compromise between the quest of a locally accurate parameterization 
of the target potential V{x) by means of a globally meaningful expansion, 
and the need of analytical expressions of its matrix elements. In particular, 
the expansion ([9]) satisfies the minimal requirement of being lower-bounded 
everywhere, provided that the last coefficient aN^^x is positive. This fea- 
ture holds independently of Nmax, in contrast to the traditional expansion in 
powers of x where the truncation to odd N^ax leads to lower-unbounded po- 
tentials regardless of the sign of the coefficients [3]. We will show in specific 
applications that the chosen form ([9]) is sufficiently flexible to approximate 
a realistic molecular potential V{x) substantially better than Vm{x), and. 
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as long as V{x) is not ill behaved, in principle arbitrarily well on a very 
wide energy range from the bottom of the well to dissociation. The Morse 
parameters and the coefficients Oj are determined by fitting the computed 
adiabatic potential V{x), thus leading to a fully ab-initio method to com- 
pute the molecular vibrational spectra, as demonstrated in a few examples 
in Sec. El Once the parameters are determined, by taking advantage of the 
algebraic properties of Vd{x), in Sec. 12.11 we build the explicit matrix repre- 
sentation of the Hamiltonian 

H = Hm + Va{x) (10) 

on a variational basis. We then perform a full diagonalization, thus obtaining 
eigenvalues and eigenfunctions for this quantum mechanical problem. In 
Sec. 12.11 we consider the basis of Morse eigenstates (jll). In Sees. 12.31 and [2^ 
we introduce an extended basis, capable to include states in the continuum, 
thus improving the variational convergency of the vibrational states close to 
molecular dissociation. 



2.1 Matrix elements on Morse basis 

The expression of the matrix elements of an operator of the form [e~"(^~^°)]*, 
such as appears in power terms [{"(x)]* = [e~"(^~^o) — 1]*, Eq. ([H]), on the Morse 
basis (jl]) are available in the literature [23j. For convenience we report here 
the analytical expression for these matrix elements computed by exploiting 
basic properties of Laguerre polynomials: 



/ 1 \n+m 

"(^-o)]t|^) = /vAr[2g _ ^ _ ^ + 1] X (11) 

(2s + Ij^ 



n — m — 1 + t\ fm — n — 1 + A /2s— 1 — m — n + ^ + t 



min(n,m) 



V n — I \ m — i 

£=0 



where Nn,Nm are the normalization constants of Eq. (|o]o. The matrix ele- 
ments of a term v{xy, are obtained by combining expressions of this kind. 
Note that the matrix element ( ITTi) is non vanishing for any n, m: therefore 
the resulting matrix is not sparse. Equation ( ITTl) permits to construct the 
matrix representation of the Hamiltonian (fTOl) . which is exact within the 
space of Morse bound states. The molecular problem is thus reduced to 
the diagonalization of a ([s] + 1) x ([s] + 1) matrix: clearly it yields exact 
eigenvalues and eigenvectors in the trivial limit H Hm (i-e. Vd{x) = 0). 



^ Equation (fTT|) holds for real t > 0, but we are interested only in integer t: in this 
case, the first non vanishing term occurs for £ = {1 — 6n.rn) max(0, min(n, m) + I — t). 
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Figure 1: The projection ^fLo l(V^jlv^A:)P of the k^^ eigenstate on the sub- 
space spanned by the Morse basis, for increasing relative importance R of the 
deformation Vd{x) added to a Morse potential characterized by [s] + 1 = 9 
bound states {h = fi = 1, s = 8.34). 

When Vd{x) ^ 0, H mixes states within the Morse basis to states orthog- 
onal to it. The matrix representation of H is therefore only approximate, 
and provides a variational extimate of the energy of the lowest bound states. 
Due to the imperfect overlap of the subspace of Morse bound states to the 
subspace of exact bound states of H, we expect that when the correction 
Vd{x) to the Morse problem becomes important in Eq. ffTOj) . the accuracy 
of the description deteriorates. This point is studied in detail in the next 
section. 

2.2 The incompleteness of the Morse basis 

Morse potential bound states are clearly incomplete for generating the full 
infinite-dimensional Hilbert space: the un-normalizable wavefunctions of the 
continuum should also be included to complete it. We expect that the in- 
complete Morse basis should provide a good description of the actual spec- 
trum for potentials weakly deformed from the Morse shape (small coefficients 
ki| ^ K))) especially for the lowest bound states. For substantial distortion 
from the Morse shape, the incomplete basis could fail badly, especially for 
states close to dissociation, and even for intermediate levels. 

The incompleteness problem is illustrated quantitatively in Fig. [1], which 
reports the total projection Yl^j=o K'^jlv'fc)^ finite-differences numerical 

eigensolutions (fkix) of the Schrodinger equation associated to the potential 
V{x) = VM{x) + a4^ [v{x)]^ to the Morse basis (jl]). As expected, the projection 
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Figure 2: The Schrodinger bound states of the anharmonic oscillator of po- 
tential energy V{x) = Vm{x) + [f (x)]"^ (sohd curve). Comparison between 
exact eigenenergies E^^ (solid levels) and eigenvalues E'^^^ obtained by diag- 
onalization of the H matrix within the Morse incomplete subspace (dashed 
levels). The good accuracy of the lowest levels deteriorates for those close 
to dissociation, and in particular the highest bound state n = 9 is obvi- 
ously missing. The pure Morse well Vm{x) is also outlined for comparison 
(dot-dashed curve). 

of the exact eigenstates of H on the subspace spanned by the Morse basis 
(jll) gradually reduces from unity (especially for the highest excited states) as 
the relative weight 




(here R = ai/Vo) of Vd{x) increases. The increase in well depth caused by 
increasing R produces an increase in the number of bound states from 9 
(pure Morse), to 12, for R = 0.5, and this makes the problem of overlaps 
even more severe. Note that the lowest eigenvectors retain a projection near 
unity even for large R. This indicates that the lowest states are determined 
fairly accurately even within the substantially incomplete basis of a strongly 
deformed Morse oscillator. In contrast, the accuracy of the states close to 
dissociation deteriorates rapidly. 

Figure [2] illustrates the effect of basis incompleteness on the approximate 
eigenenergies E'^''^ obtained by numerical diagonalization of the matrix rep- 
resentation of H, obtained by applying Eqs. and (|TT1) . As the overlap 
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analysis suggests, low levels remain accurate, while close to dissociation the 
method relying on the Morse basis gives very poor results. 

The above calculations focus on the case Vd{x) oc [t'(x)]'^. A similar 
analysis can be carried out for any term [-0(0;)]* or linear combination thereof, 
with the same definition ( |T2l) for R and qualitatively similar results. 

2.3 Beyond the Morse basis 

To improve the accuracy of the exact diagonalization method close to dis- 
sociation, it is necessary to apply it on a basis which can be systematically 
enlarged virtually to cover the entire Hilbert space. In recent years, much 
work has been done to treat at least in an approximate manner the continuum 
part of the spectrum, mostly in the context of molecules. The main reason 
is that energy regions near dissociation are becoming available experimen- 
tally: in such regions the interaction between the highest-energy bound states 
and the lowest-energy section of the continuum spectrum is supposed cru- 
cial. The quest for manageable approximations of the continuum spectrum 
has followed different routes, including analytical methods (wavefunctions 
manipulation in coordinate representation pll 125]). to essentially algebraic 
treatments related to the dynamical-symmetry viewpoint [26] or of super- 
symmetric (SUSY) quantum mechanics derivation [27j. 

An algebraic approach is the most suited for a subsequent generalization 
to potentials of the kind of Eq. (Q. Moreover, if an algebraic form of the dy- 
namical variables is available, the computation of matrix elements and other 
interesting quantities can be done by exploiting simple algebraic relations. 
We sketch here the main ideas followed in the literature, with particular 
attention to what is useful for our specific work. 

The basic idea is to build appropriate generalized ladder operators con- 
necting the different wavefunctions of a given family, and subsequently ex- 
press all relevant quantities in terms of these operators. The ladder operators 
on one side allow the identification of the dynamical symmetry of the phys- 
ical model under consideration, and on the other side lead to an algebraic 
expression for all relevant dynamical variable, including the Hamiltonian. A 
similar construction was made for the Morse potential [2H], starting from the 
Morse bound-state eigenfunctions (jlj). 

Eventually, we need a basis where suitable ladder operators act: this 
basis must be complete in C"^ [{0, 00), dy/y]. To this purpose we introduce 
the following family of wavefunctions [20] 




(13) 
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labeled by a real parameter a > 0. Following Ref. we refer to as 
quasi number state. Formally, these functions resemble Morse bound states 
(jll). However they are substantially different. In particular for any a > 0, the 
quasi number states do form a complete orthonormal set in C^[{0, oo), dy/y] 
|26j . This same set was introduced in Ref. [29j, on the basis of considerations 
related to SUSY quantum mechanics, and applied to the context of molecular 
spectra in Ref. [301 EI] • 

The construction of generalized ladder operators realized for Morse bound 
states can be repeated starting from Eq. (fT3ll [26]. The three resulting oper- 
ators are 

K_ = {ai + h) 

K+ = (ai + h) 
Ko = (tI + h, 

where hcj)^ = ncj)'^. K^, and Kq act as generalized step up, step down 
and number operator repeatedly on the states represented by Eq. (fT3ll : 

k_\rn) = c^\rn-i) (15) 

KM = C^n+llC+l) (16) 

koK) = i^ + n)\rn), (17) 

where 

Cn = ^/n{n + 2a~l). (18) 

The definition (iMl) can be inverted, to obtain the operatorial form of mo- 
mentum and Morse coordinate 

P. = '^(K-k^) (19) 

y = 2ko-{k+ + k.). (20) 

Accordingly, the whole Hamiltonian ([T]) may be expressed in algebraic form, 
as done in Ref. |2Bj: the matrix elements {(j)^\HM\4>n) "^^^ computed 
algebraically. The only non-vanishing matrix elements are those for which 
m = n, n ± 1: the Morse Hamiltonian matrix is tridiagonal on the basis 
\(pn)- All the above properties hold regardless of the value of a, for which a 
convenient choice will be made in Sec. 12. 4[ 

An instructive similarity to the algebraic structure of the harmonic oscilla- 
tor may be obtained by following the SUSY approach: in SUSY quantum me- 
chanics the Morse potential satisfies the condition of being shape-invariant, 



-2-% r^^''^"2+/^^^J 

-y+y4- (=icji + n)-y--^p^ (14) 
2 ^dy \ ^ '2 ha J ^ ' 
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so that its eigenvalues and eigenvectors can be obtained in algebraic form. 
This leads to the construction of Morse generalized ladder operators of the 
form 

%) = ^i-\ + i^p^. (21) 

parameterized by the real quantity q. The Morse Hamiltonian ([T]) can be 
expressed in factorized form 

-2 2 

Hm = ^[A\s)A{s)-s'] , (22) 



in terms of the operators of Eq. (1211) with q = s, and the formal analogy 
with harmonic oscillator algebra is evident. Define l^'g) the only state 
annihilated by A{q): A{q)\'^l) = [29] . For q = s, l^'g) coincides with the 
ground state of Morse Hamiltonian, Eq. (jlj). The other excited states can be 
obtained algebraically as 

n 

\K) = Fnl[^\s-n + i)\%--) , (23) 

£=1 

where Fn are normalization constants. Equation ( !23|) is reminescent of the 
ladder representation of the harmonic oscillator with two important differ- 
ences: (i) the excited state l^'^) is generated starting from the ground state 
of a different but related Morse problem; (ii) the ladder operators involved 
depends explicitly on the step number. 

Due to the particularly simple dependence of the operators A, on 
the parameter q, different choices of the step number q dependency generate 
new algebraic families of wavefunctions, which do not coincide with Morse 
eigenstates, but have similar algebraic properties. In particular, setting q = 
a + n, starting from the ket |(/)q) = |\E'q) annihilated by A{a), the following 
family of wavefunctions is generated: 

n 

\rn) = G:j[A\a-i+i)\ro) , (24) 

where are appropriate normalization constants. By comparing Eq. f l^T]) 
with Eq. (IT^ . the following relations emerge: 

Aia + n)rn = K_rn (25) 
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This means that, for this special (n-dependent) choice g = cr + ra, the A, A'' 
operators assume the structure of the family of ladder operators and 
the states (Elj) coincide with those of Eq. ( fT3l) . 

These SUSY considerations allow us to see the complete basis of Eq. ( fT3l) 
as a natural algebraic modification of the Morse basis, with the nontrivial 
advantage that all matrix element of operators expressed as functions of A{q) 
and A^{q) (see Eqs. (fT^ and (1201) ) can be computed by purely algebraic 
means. We will take advantage of this feature in the next section to evaluate 
the matrix elements of all operators involved in the expansion ([9]). 



2.4 Matrix elements on the quasi- number basis 



In analogy to Sec. 12.11 we derive algebraic matrix elements of powers of the 
operators e~°'^^~^°\ which appear in the expansion ([9]) of V^(x), on the basis 
fll3p of quasi number states. In terms of the operators A{q), A^{q) defined in 
Eq. fl2ip . the Morse variable y and the momentum operator can be written 
as: 



y=i2s + l)e 



-a{x—xo) 



Px = — 



2sl 

ha 
2i 



A\s) + A{s 
A{s)~A\s 



(26) 
(27) 



The A{s), A^(s) matrix elements are computed explicitly by observing 
that A{s) = A{a + n) + (s — 0" — n)I (and analogously for A"!"): 



{rjAis)m 
mA^isMi) 



Cn^rn,n-l + (s — (7 — n)5m,n 
Cn+l^m,n+l + (s — 0" — n)6m,n ■ 



i2t 



Based on the algebraic relations above, it is straightforward to obtain 
the matrix elements of all the relevant operators. In particular, the matrix 
elements of Hm are 



+ 



2Vo 
s — a 



(2s 

'n)Cm^m,,n+l + (^ 



+ (m - S + aY]5m,n 

} ■ (29) 



On the quasi number states basis (QNSB) the Hamiltonian Hm is therefore 
tridiagonal. In practice, this is a minor drawback with respect to the basis (jl]) 
of energy eigenstates, since tridiagonal matrices are diagonalized extremely 
quickly. 

We are now ready to choose the value of the a parameter: from Eq. ( |29l) 
it is evident that if o" — s equals a positive integer d, then the off-diagonal 
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matrix element {d\HM\d + 1) vanishes: a determines the sphtting of the 
Hamihonian matrix into two blocks. It is natural to set the integer d to the 
number [s] of Morse bound states, i.e. 

(7 = s - [s] , (30) 

so that < (T < 1. Under this assumption the Hilbert space Ti decomposes 
into the direct sum of a ([s] + l)-dimensional 7i~ space (bound states (jll)), 
and an infinite dimensional Ti^ space (Morse continuum spectrum) 



n = n-®n^ (31) 



In all the following we stick to the assumption 

In the same way of Eq. (1291) . the matrix elements of e~°(^~^o) g^^-g obtained 
from Eq. fl26l) : 



The matrix of e""*^^"^") is also tridiagonal. Based on Eq. (IH^ . all matrix 
elements of every term in the expansion ^ can be computed exactly. It is 
found in particular that the matrix representing [{)(x)]* is (2z + l)-band diago- 
nal|!l For the reader's convenience, explicit expressions of the nonzero matrix 
elements of the {(j)n\[v{x)Y\(j)Tn) for i = 3 to 6 are collected in Appendix [Al 

Equation (1131) represents also the basis of Morse-oscillator-like functions 
of Ref. [Sni Elj, but with an initial choice a = |[2s] + 1. The value of 
a changes as s (taken as a free parameter) is adjusted, but it is generally 
much larger than the choice of Eq. ( l30i) considered for the QNSB of the 
present work and yielding a < 1. As a consequence, the complete set of 
wavefunctions of Ref. [SUl EI] decay much faster in approaching dissociation 
(?/ — > 0). This makes convergence of the bound states much slower in that 
approach, as illustrated in Fig. [H] for the pure Morse case, where our choice 
Eq. (|30l) only requires A^^ = 9 states to reach full convergence. Similarly slow 
convergence is observed for the distorted potential of Fig. O In addition, the 
matrix elements at the left hand of Eq. (|32|) have generally nonzero value for 
any n, m, thus leading to the computational overhead of nonsparse matrices. 

In summary, the following extremely useful properties characterize the 
QNSB ( |T3l) with the choice a = s — [s]: (i) the first [s] + 1 states generate the 
same Hilbert subspace as the Morse basis of bound states; (ii) more states can 
be included simply, with a systematic improvement of the basis completeness; 



^ Like on the Morse basis, also on QNSB matrix elements of v{xY for real t may be 
obtained analytically in close form 
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Figure 3: The QNSB with a = 9, as proposed in Ref. [301 EI]: convergence 
of the bound eigenenergies of the Morse oscillator of Fig. [T] as the size Ng 
of the QNSB ( fT3|) increases. This calculation shows that the choice of a 
value of a substantially larger than our choice Eq. (130|) produces a very slow 
convergence of the eigenvalues. 

(iii) if extended to all n > 0, it becomes to a complete orthonormal basis 
of £"^[{0,00), dy/y], allowing to expand the bound solutions of an arbitrary 
1-dimensional problem, with systematically improvable accuracy; (iv) most 
importantly, it allows to express the matrix elements of all relevant operators 
in an algebraic form. In particular, the matrix relative to a term [£i(a;)]* 
is (2z + l)-band diagonal, so that the finite expansion and thus H is 
represented by an extremely sparse matrix, which is even a band matrix! We 
propose therefore the QNSB as a convenient complete orthonormal system to 
study the general problem of a molecular oscillator, in a way similar to what 
is usually done in harmonic expansion methods, but with the advantage of 
including anharmonicity from the beginning. 

3 Applications 

3.1 Convergency with basis extension 

Initially we verify that indeed the QNSB relizes basis completion as expected. 
We start considering the potential of Eq. Q, defined by 



(33) 
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Figure 4: Morse potential (dashed) and potential expansion with Vd{x) = 
^■mElsHx)]' (solid). 

for R = 0.2, as in Fig. [21 The eigenvalues E'^''^ of the if-matrix within the 
space generated by the first A^^ states in the QNSB converge rapidly to the 
exact levels E^^, those obtained by the numerical finite-difference integration, 
represented by solid lines in Fig. [2l The algebraic levels are so accurate that 
it would be impossible to tell them apart from the solid lines of Fig. [21 As a 
quantitative measure of the error of the algebraic method, we use a worst-case 
estimator 

A^. :=^max|E,»^^-i?ri, (34) 

Vq k 

which equals Ag = 2.35% for the calculation based on the Morse basis il- 
lustrated in Fig. [21 With at least iV^ = 11 states included in the QNSB, 
the bound states reach the correct number of 9 and the A^r^ diminishes 
fairly rapidly with increasing basis size: A20 = 0.051%, A30 = 0.028%, 
A40 = 0.019%. It is possible to obtain a good accuracy of this method even 
for large values of i?, namely for potentials significantly distorted from the 
Morse shape. We check that the same form ( 133|) of the potential for R = \ 
yields the correct 14 states with an excellent A30 = 0.013%. 

Good convergence against basis completion was verified also in the pres- 
ence of multiple power terms, with similar results. For example, for V{x) = 
Vm{x) + '2-OVoY^^=3[v{x)Y , strongly distorted from the Morse shape (see 
Fig. [11), the 11 bound states are also obtained with good accuracy (A30 = 
0.11%, A40 = 0.011%). 

Other tests indicate that a A^r^ within a small fraction of percent is nor- 
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mally realized by including a number of states about twice the number of 
bound states. When high powers of f (x) are included in Vd{x) with substan- 
tial weight, a larger Ng could be required for the same accuracy. Substantially 
fewer states can be included if only few low-lying levels are addressed, and 
on the contrary, Ng must be substantially increased to reach convergence 
of the topmost bound state whenever it happens to occur very close to the 
continuum, thus associated to a very extended wavefunction. 

We have therefore established that the method based on the QNSB ( IT3l) 
is reliable and efficient: given an expansion of the actual molecular potential 
in the form ([9]), this method provides accurate eigenvalues and eigenstates, 
regardless of the relative amount of deviation from the pure Morse potential. 
We come now to the investigation of the flexibility of the expansion to 
describe realistic potentials. We consider the Lennard- Jones (LJ) potential, 
and the adiabatic potential of the H2 molecule. 

3.2 The Lennard-Jones potential 

We probe the ability of the Morse expansion to describe accurately the 
standard (12, 6) LJ potential 



Like the Morse potential, Vu is a model potential, which however reproduces 
well the large-x behavior of the actual interatomic adiabatic potential of, e.g., 
a dimer of noble gas. Due to its power-law dependecies, Vlj can be considered 
as an especially "tough" potential to approximate by a series of Morse-like 
terms. 

As the functional form of the target function V{x) is assigned, it is 
straightforward to extract a few exact properties, such as well depth and 
derivatives at the minimum, which can be imposed to the approximating 
potential Vm + Vd- Exact equality of these properties can be used to fix 
a number of constraints among the Nmax + 1 parameters which determine 
Vm + Vd- The choice of the imposed exact equalities leaves some degree of 
arbitrariness to the fitting scheme. It is possible to enhance global versus 
local accuracy, by privileging properties such as the well depth, or rather a 
number of exact derivatives at the minimum. An anlysis of this sort was 
applied in a different related context (Ref. [16]). 

To start, we set Nmax = 2, i.e. Vd = and analyze the pure Morse 
potential. By imposing equal position of the minimum, equal well depth, and 
equal value of the second derivative at the minimum, we fix all three Morse 




(35) 
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Figure 5: Approximations of the LJ potential (solid): the Morse poten- 
tial (dashed), with equal position of the minimum, second derivative in the 
minimum, and well depth; the 4*'^ order (dot-dashed) and 10*'' order (dot- 
ted) Taylor expansion of Vlj- The points A and B are characterized by 
V{x) = 1% V{xo), i.e. 99% of the dissociation energy from the well bottom. 



parameters, and obtain the dashed curve drawn in Fig. O These conditions 
are the natural choice to ensure an overall agreement, which however is not 
particularly good, as the LJ and Morse potential are the prototypes of two 
quite different behaviors. As expected, the two potentials differ the most far 
from the equilibrium distance, in the approach to the dissociation limit. 

We then use the full expansion © truncated at N^ax = 4 and N^ax = 
6. We impose the following conditions: (i) same location of the minimum 
(fix value of xq), (ii) same values of second and third derivatives at the 
minimum (fix a relation between a and Vq, and a relation involving a and 
Vq and 03), (iii) same well depth (fix Vq + These 4 conditions 

reduce the number of adjustable parameters to N^ax — 3. To determine 
these residual free parameters, we minimize the RMS difference of V{x) and 
Vm{x) + Vd{x) for approximately 10^ equally spaced x-points between the 
points A, B of Fig. [5], where V{x) = 1% V{xo), i.e. 99% of the dissociation 
energy measured from the well bottom. The best fit parameters and the 
corresponding RMS discrepancy of the fit are reported in Table [H Note 
that in both calculations the fitting potential is globally well-behaved, as 
guaranteed by a-NMax > 0- The agreement between the target potential Vlj 
and the approximating expansion Vm + Vd is very good already for Nmax = 4 
(one fitted parameter only), and the RMS discrepancy improves by about 
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max ^ ■'■^ max 



Srms 0.0046 0.00043 



a 


0.10891 


0.078138 


Vo 


2.507 


4.8698 


03 


2.124 


7.6676 


04 


0.617 


5.9625 


as 




2.8206 






0.6560 



Table 1: Best fit parameters in the expansion ([9]) with Nmax = 4, 6 to the 
reference LJ potential with S = 31, A^j = 4 (so that the well depth is unity). 

one order of magnitude for N^ax = 6 (three fitted parameters). Both fitted 
potentials are visibly indistinguishable from Vlj, and have therefore not been 
drawn in Fig. [5l The fit is essentially independent of the density of the mesh 
points, except for the choice of the cut at the repulsive side. The 99% of the 
dissociation is a reasonable compromise to address the bound states. 

This accurate fit of V{x) produces extremely good accord for the cor- 
responding Schrodinger eigenvalues, as shown in Table [2l A fairly accurate 
spectrum (A40 = 0.46%) is obtained already for Nmax = 4: the highest eigen- 
value is missing but, due to its proximity to dissociation, it is indeed very 
difficult to determine. The fit with N^ax = 6 recovers the missing eigenstate, 
and improves substantially the accuracy of the eigenvalues (A40 = 0.068%), 
determined precisely by the state closest to dissociation. Note that the abso- 
lute accuracy of the eigenstates is very close to the accuracy of the underlying 
functional fit of the potential. 

For comparison with the traditional approach. Table [3] lists the exact 
LJ eigenenergies aside with (i) the approximate spectra computed for the 

relevant Taylor expansion V{x) = Vq + Ylf=2 ~" ^o)"* truncated at 

iV = 4 and 10 (drawn in Fig. [5]); (ii) the spectrum obtained by standard 
second-order perturbation theory applied to the harmonic oscillator, based on 
the fourth-order Taylor expansion; (iii) the spectrum of the pure Morse term 
drawn in Fig. [5l It is apparent that the traditional approach can only reach 
the accuracy of a few percent on few low-lying levels, but fails completely 
starting already at moderately low overtones. In contrast, even the rough 
approximation of the pure Morse term provides level estimates within few 
percent of the Vlj levels throughout the spectrum. The expanded potential 
yields a good global description of the spectrum. It would be straightforward 
to apply the present method to the recently determined adiabatic potential 
of the A12 dimer but we prefer to consider in detail a molecule where 
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n 


T T / 

LJ (exact) 


max ^ 











-0.88237 


-0.88259 


^-0.02%) 


-0.88240 


' r\ r\r\c\ V\~/ \ 

^-0.002%) 


1 


-0.67488 


-0.67592 


^-0.10%) 


-0.67494 


' r\ r\r\ V\~/ \ 

^-0.006%) 


2 


-0.50142 


-0.50387 


;-0.25%) 


-0.50146 


^-0.004%) 


3 


-0.35948 


-0.36333 


^-0.39%) 


-0.35942 


;o.oo5%) 


4 


-0.24637 


-0.25098 


' r\ A nit/ \ 

^-0.46%) 


-0.24623 


0.014%) 


r 




-O.io92 I 


-O.i6o4o 


-0.427o) 


-O.io9i6 


O.Oll/o) 


6 


-0.09514 


-0.09751 


;-0.24%) 


-0.09520 


;-o.oo6%) 


7 


-0.05078 


-0.05040 


;o.o4%) 


-0.05100 


;-0.022%) 


8 


-0.02278 


-0.01973 


;o.3o%) 


-0.02286 


;-o.oo8%) 


9 


-0.00754 


-0.00356 


;o.4o%) 


-0.00711 


0.043%) 


10 


-0.00126 


(missing) 




-0.00058 


;0.068 %) 



Table 2: Bound-state energies for the LJ potential (S = 31, Alj = 4) 
computed by numeric finite-differences solution of the Schrodinger equation 
(units of h = fi = 1), compared to the eigenvalues of the problem associated 
to the best fit Vm{x) + Vd{x), solved by diagonalization on the QNSB in- 
cluding Ng = 40 states. All reported figures are significant. In brackets the 
difference in percent of the well depth = Au/A. 



n 


LJ 


Taylor 4 


Taylor 4 


Taylor 10 


Pure Morse 




(exact) 


(exact) 


(pert, theory) 


(exact) 


(analytical) 





-0.8824 


-0.8791 (0.3%) 


-0.8842 (-0.2%) 


-0.8824 (0.001%) 


-0.8818 (0.06%) 


1 


-0.6749 


-0.6419 (3%) 


-0.6899 (-1.5%) 


-0.6742 (0.07%) 


-0.6677 (0.7%) 


2 


-0.5014 


-0.3868 (11%) 


-0.5588 (-6%) 


-0.4913 (1%) 


-0.4833 (1.8%) 


3 


-0.3595 


-0.1094 (25%) 


-0.5079 (-15%) 


-0.3065 (5%) 


-0.32862 (3%) 


4 


-0.2464 


(positive 


-0.5542 (-30%) 


-0.0989 (15%) 


-0.20369 (4%) 


5 


-0.1593 


energy 


-0.7149 (-55%) 


(positive 


-0.10850 (5%) 


6 


-0.0951 


...) 




energy 


-0.04304 (5%) 


7 


-0.0508 






...) 


-0.00732 (4%) 


8 


-0.0228 








(missing) 


9 


-0.0075 








(missing) 


10 


-0.0013 








(missing) 



Table 3: Comparison of the exact Vlj eigenstates with those of the 4 order 
traditional Taylor expansion (dot-dashed curve in Fig. [5]) obtained by exact 
integration of the differential problem, and by 2"*^ order perturbation theory. 
The eigenvalues of the pure Morse curve (dashed in Fig. [5]) are also reported 
for comparison. The differences in brackets are in percent of the well depth 
Ed. 
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Table 4: Fit parameters for model potential ©, ^max = 2,4, 12. 
PES is less close to the LJ function. 
3.3 The H2 molecule 

We acquire the accurate ab-initio determination of the H2 adiabatic potential 
of Refs. [33] . In our parameterization, we choose to fix both the well depth 
to the very highly accurate value of = 0.17447498 Ha, and the minimum 
position Xq = 1.4011 Oo p^. The other N^ax — 1 parameters, including a 
constant setting the overall energy zero, are all adjusted by minimization 
of the RMS deviation Sums from the ab-initio points provided in Ref. [53] . 
The mesh consists of all 169 available points ranging from 0.2 oq to 18 Oq 
representing the effective BO potential for H2 as computed in Ref. [33]. The 
fits include all available points. In order to enhance the convergence in the 
bound-state energy range, the points above 1% V{xo) are weigthed 1/9 of 
those in the well region. A first fit with the Morse potential term alone yields 
a rather poor approximation of {Srms = 0.035 Ha), as shown in Fig. [6l 
We then fit the expansion with Nmax = 4 and 12, obtaining 6rms = 0.006 
and 0.00003 Ha respectively. The best fit parameters are listed in Table HI 
Figure [HI reports the resulting profiles. 

Complete exclusion of the points above 1% V{xo) produces better agree- 
ment of low-order fits, but makes higher-order fits {Nmax > 8) unstable. On 
the other hand, enforcing equal accuracy to all points (no weights) makes all 
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-17340 (0.7%) 


-17614 (-0.005%) 


6 


-14612 


-14598 


-10402 (11%) 


-14282 (0.8%) 


-14599 (-0.004%) 


7 


-11830 


-11815 


-7670 (11%) 


-11491 (0.8%) 


-11815 (-) 


8 


-9287 


-9272 


-5353 (11%) 


-8979 (0.8%) 


-9271 (0.003%) 
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-6994 


-6980 


-3451 (9%) 


-6755 (0.6%) 


-6979 (0.003%) 
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-4967 


-4955 


-1965 (8%) 


-4829 (0.3%) 


-4955 (-) 
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-3222 (-0.006%) 
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Table 5: Comparison of the computed H2 bound-state energies: Ref. [33] 
includes nonadiabatic corrections; when these corrections are left out, and all 
bound states of the splined potential of Ref. [33] are computed by numerical 
solution of the differential Schrodinger problem, the energies of the second 
column are obtained. The successive columns report eigenvalues of the pure 
Morse term and of the expansion fitted to the same adiabatic potential 
of Ref. [33] . Discrepancies to the energies in the second column are reported 
in brackets: those of less than 0.001% are replaced by a slash. 
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Figure 6: Adiabatic potential for H2, from Ref. [33] (dots). Best fit Morse 
potential with a as only adjusted parameter (dashed), and fit function ([HI), 
with Nmax = 4 (dot-dashed) and N^ax = 12 (solid). Inset: detail of the 
minimum region, where the small inaccuracy of the dot-dashed curve makes 
it distinguishable from the solid curve. 

fits much less accurate. 

The model Hamiltonian is solved on a QNSB including the first A^^ = 
60 states: the computed levels are summarized in Table [5l The obtained 
eigenenergies are to be compared with those obtained by numerical integra- 
tion of the Schrodinger problem of a third-order spline going through the 
points of Ref. [33] • These reference energies (second column of Table ED dif- 
fer from the levels obtained in Ref. [33] in non including the nonadiabatic 
corrections. The agreement improves systematically as Nmax increases, with 
a worst discrepancy of 324 cm~^ for Nmax = 4, which reduces to 7 cm~^ for 
Nmax = 12. If the state closest to dissociation is not considered, all other 
levels are computed within 2.4 cm~^ of their reference values. 

In fact the N^ax = 12 spectrum could equally well be computed on a 
QNSB of Ns = 30 states only, and the same cm~^ figures are obtained. 
The reason is that the underlying Morse term is so wide and deep that it 
could host 26 bound states, so that even the pure Morse basis would be 
almost complete for the bound states of the well at hand. The Nmax = 4 
calculation generates an extra spurious bound state whose precise energy 
location converges only very slowly with Ng, while all other states require 
only Ng = 60 to remain stable to cm~^. This is due to the narrow shallow 
Morse well underlying this expansion. In such a kind of expansion, the QNSB 
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is absolutely necessary: the Morse basis would produce very poorly converged 
spectra. 

4 Discussion and conclusions 

In this work we propose a technique to compute molecular vibrational spectra 
including anharmonic effects using an appropriate virtually complete Hilbert- 
space basis, the QNSB ( |T3l) . and a smart prescription for the form of the 
model potential, namely Eq. ([9]). The substantial advantage of the proposed 
combination of Hamiltonian and basis resides in the fully algebraic form of the 
matrix elements, which can be computed quickly and efficiently, and in the 
sparse form of the final matrix. The model potential of Eq. © is sufficiently 
flexible to realize a good compromise between local and global accuracy, with 
relatively low order expansions. In principle, the virtually complete QNSB 
allows to determine accurate vibrational levels in the whole bound-states 
energy range, and even to study dissociation phenomena [^. 

In the past, a similar approach was applied to restricted forms of the ex- 
pansion ([9]). For example, the PMO prescription [Ml |l7] exclude the 
term and determine the Morse parameters a and Vq in order to satisfy the 
conditions that the second and third derivatives of the two potentials coin- 
cide at the minimum. This fixes rigidly the global shape of all other terms, 
making that expansion much less fiexible and accurate. It is eventually still 
possible to recover a good accuracy in that approximation, but it is neces- 
sary to extend the expansion ([9]) to much higher order. In fact, there is no 
special reason to fix the pure Morse contribution to some local property of 
the addressed potential, such as derivatives at the minimum, also because 
these local properties could even not be available in the realistic case of a 
microscopic potential derived by ah initio numerical calculations. 

Other algebraic methods [191 EOl EH |22] do not rely on an explicit poten- 
tial, and this makes the calculation of matrix elements of arbitrary operators 
and roto-vibrational couplings impossible to carry out without the introduc- 
tion of additional ad hoc parameters. The present method can instead gener- 
ate all matrix elements of = e~"*^, of a;*e~° ^, and of by straightforward 
algebraic manipulation. 

The method presented here is not the only one to involve smart choices 
of expansion variables. For example, in the direct-potential-fit methods 
[3Sl ESI EZl EHl EHl SDl SI] the roto-vibrational spectroscopy of diatomic 
molecules is addressed by taking the experimental spectra as input. Nona- 
diabatic effects can also be included. The validity of the potential param- 
eterization generated by such methods is usually restricted to the energy 
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range where data are available. The method of the present work generates 
instead purely vibrational spectra based on first-principle determination of 
the adiabatic PES, thus with relatively constant accuracy through the whole 
bound-state energy range. The PES parameterization proposed here is based 
on ab-initio electronic-structure calculations. It would be difficult (and of lit- 
tle interest) to attempt to use the present formalism to fit experimental lines, 
as in the direct-potential-fit approach. An extension to address nonadiabatic 
and rovibrational effects could be developed in future work. 

When this method is extended to a polyatomic context, the total basis 
size and its sparseness are the main computational issues to consider. If 
in the present scheme spectroscopic accuracy can be obtained with as few 
as 30 states per degree of freedom, the total basis size could still approach 
30^ ^ 10^ states for a tetraatomic molecule. This size could exceed the com- 
putational power available, even if the matrix to diagonalize is sparse. As 
common practice in quantum chemical calculations, a truncated basis could 
be selected by accepting slightly less accurate spectra. This truncation may 
target specific ID oscillators: in this case, when A^^ is reduced to the number 
[s] + 1 of the underlying Morse term, the QNSB and the Morse basis are 
equivalent and yield equivalent spectra. If further reduction is needed, the 
convergency of the two resulting spectra shows a characteristically different 
deterioration. Truncation of the Morse basis affects greatly the high-energy 
states near dissociation: these disappear rapidly one by one, without signifi- 
cant modification of the low-energy levels. On the contrary, the reduction of 
QNSB size affects eigenenergies throughout the whole spectrum, but with a 
much slower decrease of the number of bound states. Accordingly, the two 
basis can be both used with success, keeping in mind their different behavior 
in different situations. Alternatively, standard truncation of the full many- 
oscillator basis, based on the relative perturbative weight of the states 0, 112] 
is also made efficient by the very sparse structure of the matrix elements. 

Extension of the method presented here to vibrational spectra of poly- 
atomic molecules is not especially more complicated than the formally sim- 
ilar traditional harmonic-oscillator based expansion. In particular it lends 
itself naturally to both diagonalization (Lanczos / Davidson) [21 El El 113] 
and Green-function-based [H E] iterative methods of solution. Special care 
must be put in making sure that the total polyatomic adiabatic potential 
surface (e.g. fitted on ab initio data) has a single absolute minimum. This 
global property is easier to implement than on traditional power-series meth- 
ods where the global shape of the approximate potential is dictated by local 
properties at the equilibrium geometry. A separate treatment of "soft" tor- 
sional coordinates may be needed, exactly like within the traditional normal- 
coordinate harmonic-oscillator based schemes. Once the extension to the 
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polyatomic case is realized, the present method will compete successfully 
with the traditional methods, to yield ab-initio vibrational spectra of accept- 
able accuracy up to the high-overtones near-dissociation energy range. 



A Matrix elements on the QNSB 

We report here the nonzero matrix elements of the cubic term 

{<l>nmx)r\'Pn-s) = -^'^^^ (36) 

Mvix)m.,) = (37) 

(0„|[r)(a;)]3|0„_i) = -^{12n' - 24n + 13 + Cl_, + Cl + C^^, (38) 

+12[s]2 - 24n[s]) 

(0„|[i)(a;)]3|</.„) = (j^[C^(6n - 6^ - 5) + Cl^,{6n - 6[s] - 1) (39) 

+i2n-2[s]-lf]. 

Similarly, the nonzero matrix elements of the quartic term [^(a;)]^ read: 

mvm^n-d = ^"'(i^f;:'^" (40) 

{<PnMx)] \(pn-3) ^ ^^^^^-^4 (41) 

Mv{^)r\<Pn-2) = jf^{58 - 72n + 24n2 + Cl_, + Cl_, + Cl (42) 



{<t>n\[v{x)\'\<t>n-i) = -(i^[-10 + 26n - 24n2 + Sn^ - C^^, (43) 
+2nC2^i + Cl_^{2n - 3 - 2[s]) + 2Cl{n - 1 - [s]) - 26[s] + 48n[s] 
-24n2[s] - 2Cl^^[s] - 2A[sY + 2An[sf - %[sf] 

{ct>n\[v{x)Y\<i>n) = jT^{C^ + Ct^i + (1 - 2n + 2[s]r + Cl^, X (44) 
(2 - 8n + 24n2 + + (8 - 48n)[s] + 24^2) + Cl[Cl_^ 
+2(9 - 20n + 12n2 + Cl^^ + 20[s] - 24n[s] + 12[s]2)]}. 
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The nonzero matrix elements of the term read: 



(0,|[i)(x)]5|0„_3) = -^ii^igL^(l70- 160n + 40n2 + Ct3 (47) 
+C^_2 + C^_, + C2 + + 160[s] - 80n[s] + m[s]') 



{(^)n\[v{x)]^\(^>n-2) = -^^^[330 - 580n + SGOn^ - 80n=^ + 13C^ (48) 
-lOnCl + 9C^^i - lOnC^^i + 580[s] - 720n[s] + 2A0n'^[s] + WCl[s] 
+10Cl+^[s] + 360[s]2 - 240n[s]2 + 80[s]^ + Cl_^{17 - lOn + 10[s]) 
+Cl_,(21-10n + 10[s])] 



(0„|[i)(x)]5|0„_i) = -^j^[121 - 400n + 520n2 - 32071^ + SOn" + 

+ 14C2^, - 40nC2^i + ^On'Cl^, + C^+, + C^^^^C^^^ + 400[s] 
-1040n[s] + 960n2[s] - 320n^[s] + 40C2^Js] - 80nC^^i[s] + 520[s]2 
-960n[s]2 + 480n2[s]2 + 40C2+Js]2 + 320[s]3 - 320n[s]^ + 80[s]^ 
+2C^(21 - 40n + 20^2 + C^^^ - 40(n - l)[s] + 20[s]2) + C2_,(94 
-120n + 40^2 + C^_2 + 2Cl + C^^^ + 120[s] - 80n[s] + 40[s]2)] 

(0„|[^(x)]5|0„) = -^{C^^,(l - 10n + 10[s]) (50) 
+C^(9 - lOn + 10[s]) + C^^i[2 - 20n + 40/1^ - 80n=^ + 20(1 - 4n 
+12n2)[s] - 40(6n - l)[s]2 + 80[s]3 + Cl+^C^Ols] - 3 - lOn)] 
+C2{C2_i(13 - lOn + 10[s]) + 2[29 - 90n + lOOn^ - 40n=^ 
+10(9 - 20n + 12^2) [s] - 20(6n - 5)[s]2 + 40[s]=^ 
+C^+i(5 - lOn + 10[s])]} + (1 - 2n + 2[sm. 

The nonzero matrix elements of the term [■O(x)]^ read: 

{(pn\[v{x)\"\(t)n-6} = . (51) 
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{ct>n\[v{x)f\<Pn-.) = 12^n-4C„-^C„_.C„-,C. ^3 _ ^ ^ ^^^^ ^^2) 

(0„|[{;(x)]6|0„_4) = '^"-^^-gr^''" (395 - 300n + GOn^ + C^^ (53) 
+C2_3 + + + C2 + C^^, + 300[s] - 120n[s] + 6Q[s]') 



(0„|[^(x)]6|0„_3) = [-380 + 510ri - 240n^ + AOn^ (54) 

-6C^_, + 3nC^_, - 5Cl + 3nCl - ACl^, + 3nC^+i + C^3 x 
(-8 + 3n - 3[s]) + C^_2(-7 + 3n - 3[s]) - 510[s] + 480n[s] 
-120n2[s] - 3C^_Js] - ?,Cl[s] - 3C^+Js] - 2AQ[sf 
+120n[s]2 - 40[s]^] 



(0„|[{;(a;)]6|0„_2) = ^^^[1771 - 3960n + 3480^2 - 1440^3 + 240n4 (55) 
+CU + CU + 107^2 - IhQnCl + mn^Cl + 0^ + 590^^, - lOSnCl^, 
+60n'Cl^, + 2ClCl^, + Ci^, + Cl^.Cl^^ + 3960[s] - 6960n[s] 
+4320n2[s] - 960n3[s] + 156C2[s] - l2{)nCl[s\ + 108C,2,+Js] - 120nC2+Js] 
+3480 [s] 2 - 4320n[s]2 + 1440^2 [s] 2 + mCl[sf + eOC^^Js]^ + 1440 [s]^ 
-960n[s]3 + 240[s]^ + C2_^(179 - 204n + GOn^ + 2C2 + C^^^ + 204[s] 
-120n[s] + 60[s]2) + C2_2(275 - 252n + GOn^ + d^.^ + 2C2_, + Cl 
+C2^i + 252[s] - 120n[s] + 60[s]2)] 



{(t)n\[v{x)f\(t)n^i) = -(iq^[-91 + 363n - GOOn^ + 520n=' - 240n^ + 48n^ (56) 
-11^2^, + 42nC2+i - 60n2C2^i + AOn^Cl^, Cf,^, + 3nCt+, 
+3nC2+iC2+2 + CUi-^ + 3n - 3^) + 3C^(n - 1 - [s]) - 363[s] 
+ 1200n[s] - 1560n2[s] + 9mn^[s\ - 240n^[s] - 42C2^Js] + 120nC2^Js] 
-120n2C2+Js] - 3C^+Js] - 'iCl_^iCl^2[s\ - 600[s]2 + 1560n[s]2 - 1440n2[s]2 
+480n3[s]2 - 60C2+i[s]2 + 120nC2^Js]2 - 520[s]3 + 960n[s]3 
-480n2[s]3 - 40C2^Js]3 - 240[s]4 + 240n[s]4 - A8[sf + Cl_^ x 
[-153 + 282n - 180^2 + 40n-^ - 3C2^i + 3nC2_^i + 2Cl{-A + 3n - 3[s]) 
+3C2_2(n - 2 - [s]) - 282[s] + 360n[s] - 120n2[s] - 3C2^Js] - 180[s]2 + 120n[s]2 
-40[s]3] + 2C2(-23 + 63n - 60^2 + 20n3 + C2^i(-2 + 3n - 3[s]) 
-3(21 - 40n + 20n2)[s] + 60(n - l)[s]2 - 2Q[sf)] 



26 



= ^{Cl + Cl^, + (1 - 2n + 2[s\f + C^+i{2C^+2 (57) 
+3[1 - 4n + 20n2 + (4 - 40n)[s] + 2Q[sf\} + C^[2C2_i + 3(17 - 36n + 20/1^ + 
+(36 - 40n)[s] + 20[s]2)] + C2{179 - 696n + lOSOn^ - 800n=^ + 240n^ + C^_^ 
+3C^+i + 696[s] - 2160n[s] + 2400n2[s] - 960n=^[s] + 1080[s]2 - 2400n[s]2 
+1440n2[s]2 + 800[s]3 - 960n[s]3 + 240[s]^ + C^_i(107 - 156n + GOn^ 
+C2_2 + 2C2^i + 156[s] - 120n[s] + m[sf) + 2C2^J19 - 60n + GOn^ 
+Cn+2 - 60(2n - l)[s] + m[sf]} + C2^i{3 - 24n + 120^2 - IGOn^ + 240n^ 
+C^+2 + 24(1 - lOn + 20n2 - AQn^)[s\ + 120(1 - 4n + 12n2)[s]2 - 160(6n - l)[sY 
+240[s]4 + Cl^^i^l + 36n + GOn^ + C^^g - 12(3 + 10n)[s] + 60[s]2]}}. 
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